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ABSTRACT. Resonant motions of integrable systems subject to perturbations may 
continue to exist and to cover surfaces with parametric equations admitting a formal 
power expansion in the strength of the perturbation. Such series may be, sometimes, 
summed via suitable sum rules defining C°° functions of the perturbation strength: 
here we find sufficient conditions for the Borel summability of their sums in the case 
of two-dimensional rotation vectors with Diophantine exponent r = 1 (e.g. with ratio 
of the two independent frequencies equal to the golden mean). 



1. Introduction 

In the paradigmatic setting of KAM theory, one considers unperturbed motions <p — ► <p + ujq (T)t on 
the torus T d , d > 2, driven by a Hamiltonian H = Hq(T), where I £ Wl d are the actions conjugated 
to (p and u>o(I) = 9iHo(T). Standard analytic KAM theorem considers the perturbed Hamiltonian 
H e = Hq(T) + sf((p,T), with / analytic, and a frequency vector u)q which is Diophantine with 
constants Co and r (i.e. |u>o ■ v\ > Cq\v\~ t \/v £ Z d , v ^ 0) and which is among the frequencies 
of the unperturbed system: uj = Wo(Io) f° r some I . Suppose i?o(I) = I 2 /2 and f(f,T) = f{ip) 
for simplicity, that is assume that Hq is quadratic and the perturbation depends only on the angle 
variables. Then for e small enough the unperturbed motion ip — ► ip + oJot can be analytically 
continued into a motion of the perturbed system, in the sense that there is an e-analytic function 
h £ : T d — > T d , reducing to the identity as s — > 0, such that i/>+h s (i/>) solves the Hamilton equations 
for H ei i.e. ip + e9 v /(<p) = 0, if if) is replaced by ip + u)t, for any choice of the initial data ij). The 
function h e (called the conjugation) can be constructed as a power series in e (Lindstedt series) and 
for e small convergence can be proved, exploiting cancellations and summation methods typical 
of quantum field theory, [GBG]. We can call this the maximal KAM theorem, as it deals with 
invariant tori of maximal dimension. 

The same methods allow us to study existence of perturbed resonant quasi-periodic motions in 
quasi-integrable Hamiltonian systems. By resonant here we mean that ujq satisfies 1 < s < d — 1 
rational relations, i.e. it can be reduced to a vector (a>, 0), us £ M d s , with rationally-independent 
components, via a canonical transformation acting as a linear integer coefficients map of the angles. 
In this representation we shall write ip — (a, (3) denoting by a the "fast variables" rotating with 
angular velocity oj and by f3 the fixed unperturbed angles ("slow variables"). 

The study of resonant quasi-periodic motions is mathematically a natural extension of the max- 
imal KAM case and physically it arises in several stability problems. We mention here celestial 



mechanics, where the phenomenon of resonance locking between rotation and orbital periods of 
satellites is a simple example (as in this case the resonant torus is one- dimensional, i.e. it describes 
a periodic motion) . In general resonant motions arise in presence of small friction: the most un- 
stable motions are the maximally quasi-periodic ones (on KAM tori). In presence of friction the 
maximally quasi-periodic motions "collapse" into resonant motions with one frequency rationally 
related to the others, then on a longer time scale one more frequency gets locked to the others and 
the motion takes place on an invariant torus with dimension lower than the maximal by 2, and 
so on. Periodic motions are the least dissipative, and eventually the motion becomes maximally 
resonant, i.e. periodic. This is a scenario among others possible, and its study in particular cases 
seems to require a good understanding of the properties of the resonant motions of any dimension. 

The first mathematical result is that also in the resonant case, if u> is a (d — s)-dimensional 
Diophantine vector, a conjugation h e can be constructed by summations of the Lindstedt series; 
however the conjugation h e that one is able to construct is not analytic, but, at best, only C°° 
in e: in general it is defined only on a large measure Cantor set £ of e's, £ C [— So, £o] f° r some 
positive eo, so that C°° has to be meant in the sense of Withney. It is commonly believed that a 
conjugation which is analytic in a domain including the origin does not exist in the resonant case. 
Moreover, contrary to what happens in the maximal case, not all resonant unperturbed motions 
with a given rotation vector lj appear to survive under perturbation, but only a discrete number of 
them. This is not due to technical limitations of the method, and it has the physical meaning that 
only points /3o which are equilibria for the "effective potential" (2n) s ~ d f da f(a, f3) can remain 
in average at rest in presence of the perturbation. 

The following natural (informal) question then arises: "where do the unperturbed motions cor- 
responding to initial data (a,/3), (3 ^ /3 , disappear when we switch on the perturbation?" An 
intriguing scenario is that the tori that seem to disappear in the construction of h e actually "con- 
dense" into a continuum of highly degenerate tori near the ones corresponding to the equilibria 
/3q. It is therefore interesting to study "uniqueness", regularity and possible degeneracies of the 
perturbed tori constructed by the Lindstedt series methods (or by alternative methods, such as 
classical Newton's iteration scheme). 

In the present paper we investigate hyperbolic (sec below) perturbed tori with two-dimensional 
rotation vectors for a class of analytic quasi-integrablc Hamiltonians. Informally, our main result is 
that the hyperbolic tori, which survive to the switching of the perturbation and are described by a 
function h e that we construct explicitly, are independent on the procedure used to construct them 
iteratively (which is not obvious, due to lack of analyticity). Moreover, if r\ = ^/e, the conjugation 
is Borel summable in n for e > 0. The conjugation constructed here is the unique possible for our 
problem within the class of Borel summable functions. Of course this does not exclude existence 
of other less regular quasi-periodic motions, not even existence of other quasi-periodic solutions 
which admit the same formal power expansion as h e . 

In order to make our results more precise, we first introduce the model and summarize the results 
about existence and properties of lower-dimensional perturbed tori as we need in the following. 
Then we briefly recall the definition and some key properties of Borel summable functions, and 
finally we state more technically our main result. 

1.1. The model. Consider a Hamiltonian system 

H = ±A 2 + ±B 2 +T, 2 /(a,/3), (1.1) 
with I = (A, B) e R r x Mf the action coordinates and ip = (a, (3) E T r x T s the conjugated angle 
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coordinates. Let A = w, B = be an unperturbed resonance with \u ■ i/\ > Cq\is\~ t for v e Z r , 
v 7^ and for some C ,r > 0, and let ((A ,0), (if) + ujt,f3 )) be the corresponding unperturbed 
resonant motions with initial angles a = i/>, (3 = f3o. 
The following result holds. 

Proposition 1. [GG1] [GG2] Let (3q be a non-degenerate maximum of the function fo{(3) = 
(2ir)~ r J da /(a,/3), i.e. dpfo(flo) = and <9|./o(/3o) < 0, and let lj e M. r be a Diophantine 
vector of constants Co,t, i.e. \u> ■ v\ > Co\v\~ T for all v € HI , v ^ 0. Then for e > 0, setting 
rj = sfe, there exists a function ip —* h(ip,ri) = (a(?/>), b(i/>)) 7 vanishing as rj — > and with the 
following properties. 

(i) The functions t — > (p(t) = (a(t),/3(t)) = (ip + wt,/3o) + h(ip + ujt 7 r]) satisfy the equation of 
motion ip — —ed v f(ip) for any choice of tp G T r . 

(ii) The function h is defined for (ip,rj) e T r x [0, r/o] and can be analytically continued to an 
holomorphic function in the domain V — {rj G C : Rcr/^ 1 > r]^ 1 } . 

(Hi) (The analytic continuation ofjhis C°° in n at the origin along any path contained in T> and 
its Taylor coefficients at the origin satisfy, for suitable positive constants C and D, the bounds 
|^r9*h(^,0)| < DC k k\ T , where t is the Diophantine exponent ofu). 

Remarks. 

(1) The domain I? is a circle in C 2 , centered in r/o/2 and of radius 770/2 (hence tangent to the 
imaginary axis at the origin). 

(2) A function h was constructed in [GG1] by a perturbative expansion in power series in e = rj 2 
(Lindstedt series) and by exploiting multiscale decomposition, cancellations and summations in 
order to control convergence of the series. Eventually h is expressed as a new series which is not 
a power series in e = rj 2 and which is holomorphic in T>. The procedure leading to the convergent 
resummed series from the initial formal power series in e = rf relies on a number of arbitrary 
choices (which will be made explicit in next section) and a priori is not clear that the result is 
actually independent on such choices. Another (in principle) function h was constructed in [GG2], 
with a method which applies in more general cases (see item (4) below): we shall see that the two 
functions in fact coincide. 

(3) The function h(ip, 77) can be regarded as a function of e = i] 2 , as in [GG1] and [GG2]. As such 
the bound in item (iii) would be modified into \^d^h(tp, 0)| < DC k (2k)l T , and the analyticity 
domain in item (ii) would be as described in [GG1], Fig. 1. We also note here that the exponent 
2r + 1 in [GG1], Eq. (5.29), was not correct (without consequences as the value of the exponent 
was just quoted and not exploited), as the right one is 3t: of course for r = 1 (that is the case we 
consider in this paper) the two values coincide. 

(4) In [GG2] a similar statement was proved for /3o a non-degenerate equilibrium point of the 
function /o(/3), i.e. /3 not necessarily a maximum. If /3 is a maximum the corresponding torus is 
called hyperbolic, if /3o is a minimum it is called elliptic. In the elliptic case for e real the domain 
of definition of h on the real line is T r x £, where £ C [0, £0] is a set with open dense complement 
in [0, £0] but with as a density point in the sense of Lebesgue. The reason for stating Proposition 
1 as above is that in the present paper we shall restrict our analysis to the hyperbolic case. 

(5) If r = 1, i.e. if r = 2 and u> is quadratically irrational, than the bound on the coefficients of 
the power expansion of h in 77 at the origin is |h^ fc ^ ('*/') | < DC k k\. Hence in this case it is natural 
to ask for Borel summability of h. 

(6) Even if we consider only hyperbolic tori, in the following we shall use the method introduced 
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in [GG2], because it is more general and it is that one should look at if one tried to extend the 
analysis to the case of elliptic tori. The main difference between the forthcoming analysis and 
that of [GG2] is the use of a sharp multiscale decomposition, instead of a smooth one, as it allows 
further simplifications in the case of hyperbolic tori. We shall come back to this later. 

1.2. Borel transforms. Let F(q) be a function of 77 which is analytic in a disk centered at (2^,0) 
and radius ^ (i.e. centered on the positive real axis and tangent, at the origin, to the imaginary 
axis), and which vanishes as r\ — > as rf for some q > 1. Then one can consider the inverse Laplace 
transform of the function z — ► F(z^ 1 ), defined for p real and positive and p > po by 

If F admits a Taylor series at the origin in the form F(r/) <~ Y^k=2 FkV k then the Taylor series for 
£ _1 F at the origin is 

00 k—l 

If the series in (1.3) is convergent then the sum of the series coincides with £~ 1 F(p) for p > 
real. Of course the series expansion (1.3) provides an expression suitable for studying analytic 
continuation of C~ 1 F(p) outside the real positive axis. Note that (1.3) makes sense even in the 
case the series starts from k = 1. Then given any formal power series F(n) <~ X^feLi FkV k we shall 
define Fb(p) = J2V=i (k~i)\ as tnc B° re l transform of F(n), whenever the sum defining it is 
convergent. It is remarkable that in some cases the map F — > Fb is invertible. If this is the case 
one says that the Taylor series of F is Borel summable and we also say that the function F is Borel 
summable: this is made precise as follows. 

Definition 1. Let a function 77 — ► F(q) be such that 

(%) it is analytic in a disk centered at (j^,0) and radius j^, and admits an asymptotic Taylor 
series at the origin where it vanishes, 

(ii) its Taylor series at the origin admits a Borel transform Fb(j>) which is analytic for p in a 
neighborhood of the positive real axis, and on the positive real axis grows at most exponentially as 
p — > +00, 

(in) it can be expressed, for 77 > small enough, as 

F(r])= f + e- p ^F B (p)dp. (1.4) 
Jo 

Then we call F Borel summable (at the origin). 
Remarks. 

(1) If F is Borel summable, then one says that F is equal to the Borel sum of its own Taylor 
series. 

(2) For instance one checks that a function F holomorphic at the origin and vanishing at the origin 
is Borel summable; its Borel transform is entire. 

(3) The function F(n) = X^feLi ^ T+k~ 1S n0 ^ ana hyti c a t the origin but it is Borel summable. 

(4) If F, G are Borel summable then also FG is Borel summable and (FG)b(p) = J Q P Fb(p')Gb(p- 
p')dp' = (Fb * Gb)(p) on the common analyticity domain of Fb and Gb, with the integral which 
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can be computed along any path from to p in the common analyticity domain. Of course by 
definition \Fg * Gb\ < \Fb\ * \Gb\, where in the r.h.s. the convolution is along any path from 
to p in the common analyticity domain of Fb and Gb (note however that now the convolution 
\Fb \ * \Gb\ depends on the path). 

(5) The Borel transform of rf is p k ^ 1 /(k — 1)!. The Borel transform of 77/(1 — an) is e ap . If 
F B (p) = e a Wfci! and G B (p) = e ap p k2 /k 2 l, then F B * G B (p) = e ap p kl+k2+1 /(h + h + 1)! 

(6) More generally if F i} i = 1, 2, have Borel transforms FiB bounded in a sector around the real 
axis and centered at the origin by < (7e ai l p l + ' 3i l Imp l \p\ ki /ki\, with en and Pi real, then 
I -Fib * F 2B (p)\ < C 2 e a \ p \ +fi \ lmp \ |p| fcl+fc2+1 /(fc 1 + k 2 + 1)! where a = max a, and f3 = max/3,. We 
shall make use of this bound repeatedly below. 

1.3. Main results. We arc now ready to state more precisely our main results. 

Proposition 2. Let us consider Hamiltonian (1.1) with r — 2 and f (a, (3) an analytic function of 
its arguments. Let /3q G T s satisfy the assumptions of Proposition 1 and u> G M 2 be a Diophantine 
vector with constants Co > and r = 1, i.e. \u> ■ i/\ > C 1 1>* | — 1 for all v G U ' , v 7^ 0. Then there 
exists a unique Borel summable function \\{xp,rj) satisfying properties (i)-(iii) of Proposition 1. 

Note that, once existence of a Borel summable function h(ip, n) satisfying properties (i)-(iii) of 
Proposition 1 is obtained, the uniqueness in the class of Borel summable functions is obvious, by 
the very definition of Borel summability. In fact all such functions have a Borel transform he (■»/>, p) 
that is p-analytic in an open domain enclosing R + (hence also in a neighborood of the origin), 
where they coincide (because they all have the same expansion at the origin), then they all coincide 
everywhere. 

The proof of Proposition 2 will proceed by showing Borel summability of (one of) the function(s) 
h(i/>,fy) constructed in [GG2]. In particular a corollary of the proof is that h(il>,r]) constructed in 
[GG2] is independent of the arbitrary choices mentioned in Remark (2) after Proposition 1. 

Our proof of Borel summability of h does not use Nevanlinna's theorem, [Ne] [So]. In fact we 
failed in checking that h satisfies the hypothesis of the theorem. Our strategy goes as follows. 
We introduce a sequence of approximants to h, that is naturally induced by the multiscale 
construction of [GG2]. We explicitly check that is Borel summable and that its Borel transform 
is entire. Then we show inductively that the analyticity domain of h^-* is a neighborood B of R + 
(not shrinking to as N — > 00) and that grows very fast at infinity in B (in general faster 
than exponential). However the results of Proposition 1 imply that the growth of on the 

positive real line is uniformly bounded by an exponential. Then Borel summability of h follows by 
performing the limit N — > 00 and using uniform bounds that we shall derive on the approximants 
and on their Borel transforms. 

In the next section we will recall the structure and the properties of the resummed series obtained 
in [GG2], defining the function h(i/>, 77) of Proposition 1. In Section 3 we define the sequence 
of approximants and we show that the inverse Laplace transform of is uniformly bounded by 
an exponential on the positive real line. In Section 4 we prove Borel summability of h(xp,r]) in 
the easier case in which the perturbation /(a,/3) in (1.1) is a trigonometric polynomial in a. In 
Appendix Al we discuss how to extend the method to cover the general analytic case. Finally, in 
Appendix A2 wc show that the same result applies to the function h constructed in [GG1]: this 
allows us to identify the functions constructed with the two methods of [GG1] and [GG2], since 
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they are both Borel summable and admit the same formal expansion at the origin. 

2. Lindstedt series 

Denote by a(ip),b(ip) the at, [3 components of h, respectively. In [GG2] an algorithm is described 
to construct order by order in rf the solution to the homologic equation 



f (« • ^) 2 a(</>) = V 2 d a f(iP + a(t/>), 0o + bW), 
\ ( U ■ ^) 2 b(^) = rfdpfW + a(^),/3 + b(t/>)). 



The resulting series, called the "Lindstedt series", is widely believed to be divergent. A summation 
procedure has been found which collects its terms into families until a convergent series is obtained. 
The resummed series (no longer a power series) can be described in terms of suitably decorated 
tree graphs, i.e. h can be expressed as a sum of values of tree graphs: 



0ee„, 7 



where h„ is the u-th. coefficient in the Fourier series for h, and 7 = {1, . . . ,2 + s} labels the 
component of the vector h u (recall that 2+s is the number of degrees of freedom of our Hamiltonian, 
2 being the number of "fast variables" a and s being the number of "slow variables" (3). Q v . 7 is 
the set of decorated trees contributing to h 1:V (r]) and, given 9 £ 9„ i7 , Val(0) is its value, both still 
to be defined. 

We now describe the rules to construct the tree graphs and to compute their value. We shall 
need the explicit structure in the proof of Borel summability in next sections, and this is why 
we are reviewing it here. Given the rules below one can formally check that the sum (2.2) is a 
solution to the Hamilton equations, see [GG2]. A few differences (in fact simplifications) arise 
here with respect to [GG2], and we provide some details with the aim of making the discussion 
self-consistent. Essentially, the changes consist of: (i) shifting the order of factors in products 
appearing in the definition of the values Val(0) to an order that makes it easier to organize the 
recursive evaluation of several Borel transforms; see remarks following (2.9), and (ii) using a sharp 
multiscale decomposition; see item (f) below. 

Consider a tree graph (or simply tree) with k nodes v\, . . . ,Vk and one root r, which is not 
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considered a node; the tree lines are oriented towards the root (see Fig.l). 




Figure 1. A tree 9 with 12 nodes; one has p VQ = 2,p vi = 2,p V2 = 3,p V3 = 2,p Vi = 
2. The length of the lines should be the same but it is drawn of arbitrary size. The 
separated line illustrates the way to think of the label r) = (7', 7). 

The line entering the root is called the root line. We denote by V(9) and A(0) the set of nodes 
and the set of lines in 9, respectively. 

(a) On each node v a label v v e Z 2 , called the mode label, is appended. 

(b) To each line £ a pair of labels 77 = (7', 7) is attached. 7' and 7 are called the left or right 
component labels, respectively: 7' € (1, . . . ,2 + s) is associated with the left endpoint of £ and 
76 (1, . . . , 2 + s) with the right endpoint (in the orientation toward the root, see Fig.l). The label 
7' associated with the root line will be denoted by 7(0). 

(c) Each node v will have p v > entering lines £\ , . . . , £ Pv . Hence with the node v we can associate 
the left component labels j' vl , . . . , *y' of the entering lines, and an extra label j v o = Je attached to 
the right endpoint of the line exiting from v. Thus a tensor d~, v0 j vl ...j fv v {0o) can be associated 
with each node v, with <9 7 denoting the derivative with respect to /3 7 if 7 > 2 and multiplication 
by iv 1 if 7 < 2. 

(d) A momentum vt is associated with each line I — v'v oriented from v to v': this is a vector in 
Z r defined as vt — J2 w <v Vw - ^ ne roo * momentum, that is the momentum through the root line, 
will be denoted by u{0). 

(e) A number label kg e {1, . . . , |A(0)|} is associated with each line I, with U£ e \^{ke} = {1, . . . , 
|A(0)|}. The number label is used for combinatorial purposes: two trees differing only because of 
the number labels are still considered distinct. 

(f) Each line I also carries a scale label ng — —1,0, 1,2,.. .: this is a number which determines 
the size of the small divisor oj ■ i>i, in terms of an exponentially decreasing sequence {7 P }^L °f 
positive numbers that we shall introduce in a moment. If = then = — 1. If |u> • vg\ > 70 
then ng = 0, and we say that the line i (or else u> ■ vg) is on scale 0. If 7 P < \lj ■ u\ < 7 p _i for 
some p then ng = p, and we say that the line £ (or else u> ■ vg) is on scale n. The sequence {7 P }^L 
is such that j p G Co[2~ p ~ 2 , 2~ p_1 ) for all p > and, furthermore, u> ■ v not only stays bounded 
below by C 1 1^ | — 1 (because of the Diophantinc condition) but it stays also "far" from the values 
7 P for v not too large, i.e. for \v\ at most of order 2 P ; cf. [GG] for a proof of the existence of the 



7//e6braio/2008; 15:14 



7 



sequence (without further assumptions on uj). Precisely, 

(1) \u)-v\> CoM"\ 

(2) min | \u ■ u\ - 7 P | > C 2~™ 



0<p<n 



0/i/eZ r , 

if n > 0, < \u\ < 2 ( "~ 3) . 



(2.3) 



Note that the definition of scale of a line depends on the arbitrary choice of the sequence {j p }: 
we could as well used a sequence scaling as -f p ~ -f~ p ', with 7 any number > 1 instead of 7 = 2; 
or we could have used a smooth cutoff function (as in [GG2]) replacing the sharp cutoff function 
1(7p < \uj ■ f \ < 7 P -i) implied in the definition above. 

(g) The scale labels allow us to define hierarchically ordered clusters. A cluster T of scale n is a 
maximal connected set of lines I on scale nt, with ne < n, containing at least one line on scale n. 
The lines which are connected to a line of T but do not belong to T are called the external lines 
of T: according to their orientations, one of them will be called the exiting line of T, while all the 
others will be the entering lines of T. All the external lines i are on scales n£ with ng > n. The 
set of lines of T, called the internal lines of T, will be denoted by A(T) and the set of nodes of T 
will be denoted by V(T). 

(h) Not all arrangements of the labels are permitted. The "allowed trees" will have no nodes with 
momentum and with only one entering line and the exiting line also carrying momentum. We 
also discard trees which contain clusters with only one entering line and one exiting line with equal 
momentum and with no line with momentum on the path joining the entering and exiting lines 
("self-energy" clusters or "resonances"). 

Remark. One can verify that chains of self-energy clusters can actually appear in the initial 
formal Lindstedt series. One of the main points in [GG1] and [GG2] is to show that if one modifies 
the series by descarding all chains of self-energy clusters, then the resulting series is convergent 
(a form of Bryuno's lemma that appears in KAM theory). In both [GG1] and [GG2] it is shown 
that, in order to deal with chains of self-energy diagrams one can iteratively resum them into the 
propagators (i.e. the factors associated with the tree lines in the value of a tree, see below for a 
definition) , that will then turn out to be different from those appearing in the naive formal Lindstedt 
series (which are simply (u> ■ v)~ 2 ). Such resummation is the analogue of Dyson's equations in 
quantum field theory and the iteratively modifed propagator has been, therefore, called the dressed 
propagator. Here there is further freedom in the choice of the self-energy clusters. The idea is 
that the self-energy clusters must include the "diverging contributions" affecting the initial formal 
power series. But if we change the definition of self-energy clusters by adding to the class a 
new class of non-diverging clusters, the construction can be shown to go through as well. We 
find convenient the specific choice above but this is of course not necessary. This is the second 
arbitrary choice we do in the iterative construction of the resummed series. It can fuel doubts 
about the uniqueness of the result which can only be dismissed by further arguments (like the 
Borel summability that we are proving). 

The set of all allowed trees with labels 7(6*) = 7 and i>{6) = u is denoted by 0„ 7 (this is the set 
appearing in (2.2)). The labels described above are used to define the value Val(#) of a (decorated) 
tree 9 € 0^ 7 : this is a number obtained by multiplying the following factors: 

(1) a factor F v — 9 7tj07 ' i ... 7 ^ /i/„(A)), called the node factor, per each node v; 

(2) a factor gf 1 ^ = ffy^(u; • v>e; n), called the propagator, per each line of scale m, momentum vg 
and component labels 7^,7^, see items (I)-(V) below for a definition. 
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The value is then denned as 



Val(0) = l( J] KM II ffM- ( 2 - 4 ) 

»6V(«) £eA(e) 

where it should be noted that all labels 7 (of the tensors F v and of the matrices </j™^) appear 
repeated twice because they appear in the propagators as well as in the tensors associated with 
the nodes, with the exception of the label 7 associated with the left endpoint of the line ending in 
the root (as the root is not a node and therefore there is no tensor associated with it). 

Adopting the convention of summation over repeated component labels Val(0) depends on the 
root label j(8) = 7 = 1, . . . , 2 + s so that it defines a vector in C 2+s . 

The recursive definition of the propagators is such that the series in (2.2) is convergent and gives 
the v-th Fourier component of the function h(ip,r]) in Sect. 1. The definition of propagators 
we adopt here is the same introduced in [GG2]: the definition in [GG1] is slightly different (see 
Appendix A2), but it has the drawback that it is specific for hyperbolic resonances, while the 
definition in [GG2] can be (expected to be) extended also to the theory of elliptic resonances and, 
therefore, might turn out to be useful in view of possible extensions of the main results of this 
work to elliptic resonances. 

(I) For n = — 1 the propagator of the line I is defined as the block matrix 

(0 (-^/oV))- 1 ) ' (2 ' 5) 

def 

(II) For n — 0, if the line £ carries a momentum v and if x = u? • i/, the propagator is the matrix 

g? ] =gM(x;r,)= 2 f 2M , (2.6) 

X Z + T] Z M 

with M d = ( ® 2 ^ . ] . By the assumptions of Proposition 2 one has M > 0. 

(III) For n > the propagator is the matrix 

f+Jn^y (2 ' 7) 

with M [ - n] (x;r)) = M [0] (x; 77) + M [1] (x: if) + . . . + M [n] (x; 77), where M [0] (x; 77) = r7 2 M , whereas 
MW(x;r]), j > 1, are matrices, called self-energy matrices, whose expansion in r\ starts at order 
rf and are defined as described in the next two items. 

(IV) Let T be a self-energy cluster on scale n (see item (h) above) and let us define the matrix 1 
Vt(u> • v; if) as 

v ^-^ = -jmu{ n *0( n 9 [ r ] ), (2.8) 

1 v vev(T) ieA(T) 



1 This is a matrix because the self-energy cluster inherits the labels 7, 7' attached to the left of the entering line 
and to the right of the exiting line. 
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where, necessarily, ng < n for all I e A(T). The matrix (2.8) will be called the self-energy value 
of T. The set of the self-energy clusters with value proportional to r] 2k , hence with k — 1 internal 
lines with ni > 0, and with maximum scale label n will be denoted S^ n . 

(V) The self-energy matrices A4^(x; vf), n > 1, are defined recursively for |x| < 7„_i (i.e. for x 
on scale > n) as 

oo 

^ H (^) = E E m^), ( 2 - 9 ) 

where the self-energy values are evaluated by means of the propagators on scales p, with p = 
-l,0,l,...,n-l. 

Remarks. 

(1) With respect to [GG2] the second argument of the propagators (and of the self-energy values 
and matrices) has been denoted r\ instead of e; we recall that the variable rf appearing here is the 
same as the variable e appearing in [GG1] and [GG2]. We make this choice because it is natural 
to study Borel summability in w and not in e = rj 2 . 

(2) The association of the factors rj 2 with the lines themselves rather than with the nodes (as in 
[GG1] and [GG2]) will be more convenient when considering the Borel transforms of the involved 
quantities. 

(3) The multiscale decomposition used in [GG2] may look quite different from the one we are using 
here, but this is not really so. First, even though the decomposition in [GG2] was based on the 
propagator divisors A^(x;e) = min^ \x 2 — A^"'(e)|, where the self-energies A^(e) were defined 
recursively in terms of the self-energy matrices, in the case of hyperbolic tori one has identically 
AW(i;£) = x 2 : indeed all self-energies which are non-zero are strictly negative. Then the only 
real difference is that here we are using a sharp decomposition instead of a smooth one, but the 
latter is not a relevant difference. In fact the choice of the sequence {7 P }^L implies that the 
lines appearing in the groups of graphs that will be collected together to exhibit the necessary 
cancellations have currents v such that v ■ u) stays relatively far from the extremes of the intervals 
[7p+i,7 P ] that define the scale labels, and this allows us to use a sharp multiscale decompositions 
instead of the smooth one used in [GG2] . In other words this change with respect to [GG2] is done 
only to avoid introducing partitions of unity by smooth functions and the related discussions. 

Therefore the expression (2.2) makes sense and in fact the function h mentioned in Proposition 
1 is exactly the Fourier sum of the r.h.s. of (2.2). In particular in [GG2] it was proved that the 
Fourier sum of the r.h.s. of (2.2) satisfies the properties (i)-(iii) in Proposition 1. 

3. Integral representation of the resummed Lindstedt series 

Given the definitions of Sect. 2 consider the function defined in the same way as h but 

restricting the sum in (2.2) to the trees containing only lines of scale n < N. 

The functions h<- N > have the "same" convergence and analyticity properties of the functions h 
and the same bounds on the Taylor coefficients at the origin. Moreover h: this is 

a consequence of the intermediate steps in the proof of the above proposition in [GG2], as the 
strategy of the proof is to define making sure that the convergence and analyticity properties 
are uniform in N. In fact h^ N ^ is even analytic in rj near the origin for \rj\ < t/n (but r/N j^^^' 0). 

Therefore the functions h' '• ) are trivially Borel summablc, and have an entire Borel transform, 
but the growth at p — > +oo of their Borel transforms is iV-dependent while, to show Borel summa- 



7//ei>i>raio/2008; 15:14 



10 



bility of h, uniform estimates are needed. This section is devoted to a first attempt at such bounds 
which uses minimally the informations on the resummed series that can be gathered from [GG2] , 
i.e. the convergence properties just mentioned. 

The Borel transform of the functions h( N \ij>,r]) is an entire function that can be written for p 
real and positive as 

(h^) B (^,p) = C- 1 h(^,p)= e ^hW(V,-)— , pe» + , (3.1) 

where t]n is the convergence radius of h.( N \ip,ri). The key remark is that we also know that 
by property (ii) in Proposition 1 (actually by the same property for that follows from the 
construction in [GG2]) the function h^ N \ip, -) is analytic for \z\ > 2r7 ~ 1 , so that the integral in 
(3.1) can be shifted to a contour on the vertical line with abscissa p > 2ri ( ^ 1 , i.e. with TV-independent 
abscissa. Therefore 

(hW) s (<M=/ e^hW^,-)^, (3.2) 

Jp-ioc z Z7T " 1 

and, for all N, the function h( N \ip, i) is uniformly bounded by O(rjp) on the integration contour, 
because, by property (iii) in Proposition 1, h is twice differentiable at the origin along any path 
contained in T> (in particular along the circular path Re?? -1 = p). Hence the latter boundedness 
property of h^ N \tp, j) and (3.2) imply the bound, for p > and for a suitable constant C, 

|(hW) B (^,p)| < Ce* VpeM+, (3.3) 

and for all p > 2?7 ~ 1 . The existence of the limit limAr^oo h.( N ^(ip, i) = h(ip, j) for A small (by 

[GG2]) implies existence of the limit F(ip,p) as N oo of (hW) s {tp,p) for peR + and F(tp,p) 
satisfies the bound (3.3) on M + . 
Hence the functions h(ip, 77) can be expressed, for < 77 < r?o, as 

/•OO /'OC 

h(V>,77)=/ e- p ^( lim (hW) B (ip,p))dp= / e-P/"F(-0,p)dp, (3.4) 
Jo Ar ^°° Jo 

which provides us with an integral representation of the resummed series and shows that the 
resummation (2.2) generates a Borel sum of the formal Lindstedt series provided F(ip,p) can be 
shown to be analytic in a neighborhood of the positive axis, as required by the very definition of 
Borel summability, see property (ii) in Definition 1 of Section 1.2. 

4. Borel summability 

By the remark at the end of Sect. 3 Borel summability of h will be established once the natural 
candidate for its Borel transform, namely the function F in (3.4), is shown to be analytic in a 
region containing the positive real axis. This will be done here assuming, for simplicity, that the 
perturbation / is a trigonometric polynomial. The general case of an analytic / is slightly more 
involved and will be treated in Appendix Al. 
Assume, inductively, that Vx = u ■ v the functions (M^-^)b(x',p) are entire functions of p and 

\\(M [ ^ k] ) B (x;p)\\ < r\p\e d "M, for all x on scale k' , with k' > k > 0, (4.1) 
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for all p complex. The matrix norm in (4.1) is ||.M|| = maxj ^ |-A4y|. Note that (M^)b(p) = 
pMo (so that the inductive assumption in (4.1) is valid at the first step k — with do = if 
r > ||Mo||) and (g^)B(x;p) is entire and it is given by 



{g n )B[x . p) = ^^0^ => \\{ 9 ^) B {x-,p)\\<f 2 e^, vn> 70 , (4.2) 



with c = vl | Mo 1 1/70. 

Supposing the inductive assumption to be valid for k = 1, . . . , N — lwe remark that this implies 
a bound on (g^)B(x;p) for k = 1, . . . , N — 1 via the expansion 



(.").(«) - (S E (=^)X <«> 

m=0 



Taking the Borel transform and performing all convolutions along a straight line from to p, for 
x on scale 1 < k < N — 1, we get 

m=0 m=0 v ' j (4.4) 

< N dfcbl+rV^-i^i = |p| | p | 

— 9 — 9 

with /* m d =/ *...*/ m times and c fe = d k + T 1 / 2 ^ 1 . 
Then (M [ ^ N] - M [a] ) B (x;p) is estimated via (2.9), that is 

\\(M^-MM) b (x;p)\\ < 



< 



00 1 1 

e e mm it iN^ lp| )( n ii^- 1] h)( n 

fc = 2 TPU JV_1 S' R €EA(T) « (£A(T) t>GV(T) 

C 3=0 "k,j n e >0 n e = -l 



(4.5) 



where X£ = lj ■ v>£, the f\* is a convolution product and 

\\ F v\\= , maX , EK F ")7.o,7i 1 7i a ...7i p J- ( 4 - 6 ) 
7„ii7 V 2'"-i7vp v 7tjQ 

Hence bounding |p| * (j]* He Cn * |p| ) by e^-'W |p| 2 ' £ - 1 /(2/c - 1)! and using the estimates in [GG2] 
to control the sum over the self-energy clusters, the bound becomes 

00 I |2fc— 1 

\\(M^) B (x;p)\\ <rip| + ^r 2fc e ™-iH < rble^H, (4.7) 

fe=2 ^ 

where T is a suitable constant derived in [GG1] and djv = cat_i + T. 

Remark. The step leading from (4.5) to (4.7) is non-trivial and the possibility of bounding the 
small divisors xj and the sum over self-energy clusters, after defining the propagators as above, is 
the main technical aspect of the work in [GG2] . We take the existence of T from [GG2] . We do 
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not repeat here the analysis performed in Sects. 5 and 6 (and the corresponding Appendices) of 
[GG2]: the constant T 2 has been called s^ 1 in Theorem 1 of [GG2]. 

Thus the inductive assumption holds for all N, the constants cjv, djv can be taken c2 N for some 
c, and for all x on scale N one has 

||( 5 [Arl )B(x; P )||<M e2 "cbl (4 .8) 

This leads to a bound on (h^) B (ip,p), via (2.4) and (2.2): 

00 1 1 

i(h%(^)i<EEE txL (n*>i^ w )( n n*ni)( n 

fe=i i/ e„ , 1 Wl ' «eA(») < <ea(t) t>ev(T) 

n £ >0 n f = -l 

00 I \2k—l 00 I \2k—l 

< p2fc 1PI c2»-» |p| < V- V 2k IP] c'fc |p| < r 2 1 i r|p|e<= W 

"fe (2*-l)! C (2fc-l)! 6 " lPl 

(4-9) 

because the maximum scale n max of the lines of a graph with /c lines can be at most log 2 (&fc) by our 
assumption that / is a trigonometric polynomial: in fact the maximum momentum on a line can 
be \u\ < bok (for / a trigonometric polynomial), so that the smallest x can be ^ if Co and t = 1 
are the Diophantine constants, hence the scale of x can be at most log 2 (6ofc/4) and c2" max < c'k 
for a suitable c'. 

Therefore the functions (h^ N ^)B{ip;p) are entire and bounded by 

\(h^) B (4,;p)\<T 2 \p\e r ^' lPl (4.10) 

independently of N. Since for \p\ small the functions (h^ N ^)B(ip;p) converge to F(p), simply 
because of the bound \-gd k h(xp, 0)| < DC k k\ on the Taylor coefficients at the origin of h(ip, rj), it 
follows that F(ip;p) is entire inp and everywhere the limit of (h( JV ') B (t/';p). Therefore (by Vitali's 
theorem) F(ip;p) is not only holomorphic in p near (which is a result on the Lindstedt series at 
the origin) but it is entire. In particular it is holomorphic around the real axis and therefore by 
(3.4) h(ip;r]) is Borel summable. 

5. Concluding remarks 

(1) It is interesting to stress that the recursive bounds on the Borel transform of the propagators 
in (4.4) have been derived without making use of the cancellations that played such an essential 
role in the theory in [GG2] (and [GG1]) and by "undoing" at each step the resummations which led 
to the construction of h and to the proof of Proposition 1, see (4.3) above. However the properties 
of h and the result of Proposition 1 have played a key role in the derivation of (3.3) and (3.4). 
Without the uniform bounds (in N) on the convergence radii in p of the series expressing (hW ) b , 
which depend on the cancellations and on the resummations, the bounds in Sect. 4 or, in the 
non-trigonometric case, of Appendix Al, would remain the same but they would be useless for our 
purposes of establishing (3.4) and the Borel summability. 

(2) The assumption that / is a trigonometric polynomial has been heavily exploited and an extra 
idea is necessary to deal with the more general case of analytic /: this is discussed in Appendix 
Al. 
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(3) It has been remarked above that the resummation procedure is based on several arbitrary 
a priori choices which therefore may lead to the existence of several solutions of h with the same 
asymptotic series at e = 0. All the choices in [GG2], as well as that used in [GG1] (see Appendix 
A2), lead however to a Borel summable series: this proves that all solutions coincide and the results 
of the resummations are independent of the particular choices provided the Diophantine constant 
r is t = 1 (hence r = 2 and u> is a Diophantine vector with r = 1, e.g. W1/1V2 is a quadratic 
irrational). 

(4) If r > 1, hence if r > 2, the problem of Borel summability and of independence of the result 
from the summation method remains open. 

(5) The existence or non-existence of solutions h which are C°° at the origin and give solutions 
to the equations of motion but which are not Borel summable is also an open problem. Note that 
even in the case of non-resonant motions the problem of the uniqueness at fixed lj is not trivial, 
and also the recent results in [BT] do not exclude the possibility of other quasi-periodic motions 
besides those constructed through the KAM algorithm. 

(6) The case e < with (3q a minimum point for /o(/3), i.e. of elliptic motions is quite different. 
We have used on purpose the resummation technique of [GG2] which works for hyperbolic as well 
as for elliptic resonances, but the present results still only apply to the hyperbolic case and it is not 
clear whether the above techniques can be extended to prove Borel summability of the parametric 
equations for h in elliptic cases. 

Appendix Al. Analytic, non-polynomial perturbations 

In this Appendix we want to prove real analyticity of ¥{tp,p) in the case / is a generic analytic 
function (rather than a trigonometric polynomial); in this case the function F p) will turn out 
not be entire but only analytic in a strip of width 2k > around the real axis. 

The proof of this claim will be based on an inductive assumption on (g^)B(x;p) formulated by 
introducing the matrix g^(x;n) = X 2 +-M f<„) (a ,.^ ) for all \x\ < 7„_i. Note that if Xn (x) = 1 (j n < 
\ x \ < 7n-i) is the indicator of the scale of x, the propagator g^ n \x;ri) is given by g^(x;ri) = 
Xn{x)g^{x;ri). Furthermore the matrices g^(x;r)) satisfy the recursive equations 

(g^ix^))- 1 = (g^- 1 Hx;r ] )y 1 +r 1 - 2 M^(x;n), V|x| < 7 „-i • (Al.l) 



We suppose, inductively, that for n — ^/||M || one has 

\0 n] ) B (x;p)\\ < XoMe^+^kr^^bl+^olimpll.r 1 5 V|z|<7n-i, n>0. (A1.2) 

Note that for n = the explicit expression of (g^)B given by (4.2) implies that (A1.2) is valid 
with Kq = yfs (where s is the dimension of the non trivial block in Mo) and cq — c' — 0. The 
constant Kq comes from our choice of the matrix norm ||M|| = maxj ^ i |Mj|: with this choice 
for any d x d matrix we have <i -1 / 2 ||M||2 < \\M\\ < d}/ 2 \\M\\2 where || • H2 is the spectral norm, so 
that in particular || sin M/M \ |, || cosM|| < d 1/2 . 

Assuming (A1.2) for n < N — 1 and performing all convolutions along the straight line from the 
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origin to p, we get the following bound on {M\ n], )b{x;p): 
\\( V - 2 MW) B (x;p)\\< 

<E E utW( n^oM e ^^< /2 )H-o|i mP |<) 



|A(T)|! 

n,>>0 



( n h^- 1] ii)( n ii^n) (AL3) 



<6A(T) t>ey(T) 



oo | 1 2 fc — 3 

< ^ p2fc ^ gCcjV-i+cW-iTw'-i^lPl+^ol^Phw-i e" 2 ^" < £, o | p | £ d N \p\ e -K2" 



k=2 



where Y and n are deducible from Appendix A3 of [GG2], Dq = T 3 , = cjv_i + c' N _ 1 "f^ 1 ^ 2 and, 
using that 7^_x < 4C ~ 1 2 JV , we chose |Imp| so small that 4K C ( J~ 1 |Imp| < n. 
Using (A1.3) and (Al.l) we can get a bound on ^ n ^)b{x;p)'- 

\P\ 



(9n) b (x; P )\\ < (KoMe^-i+^-xkr^^H+^limpll.r^^ 



oo , , , (-41-4) 

* E [(^o^e^^+^-^^^'IPl+^I^P^r 1 ) * (D \p\e d ^e-* 2N )} , 

k=0 X 

I |4fc+l 

and, since \p\ * (\p\ * \p\)* — for k > 0, the fc-th term in the sum is bounded by 

(K D e- K2N ) k \p\ \p\ 4k (rfw+cW _ i | x |-i/2 )b | +K0 |i mp | kr i 
(4fc + l)! ^ z 2fe ' 1 ° J 

Summing (Al.5) over A: > and comparing the result with the inductive assumption (A1.2), we 
see that we can take cat = = cjv-i + c'n-i'Jn^-i and c' N — c' N _ 1 + KoDoe~ K2 . Solving the 
iterative equations for cn,c' n , we see that, for x on scale n, (<7^)b can be bounded as 

\\(9 [n] ) B (x;p)\\ < A- M e c2«/=b| +Kl |i mp |2« ; 7n < |x| < 7n _ i ^ n > ^ (AL6) 

for some constants c, Ki . 

Plugging this bound into the expansion for h^(i/>,p), denoting by N(9) the maximal scale in 9 
and choosing |Imp| small enough, we finally get 

|(hW) B (^,p)| < 

^EEE E iot( IT §we— h-' 1 -'-) 

fc=l n >0 eee^,7 1 v *€A(o) ^ 

N(e)=n n^>0 

( n n ii^ii)^ (AL7) 

<eA(T) uey(T) 



»=-i 



oo I |2/c — 1 

<^ ^ r 2fe_lPl__ e c2"0/2|p| + Kl | Im p|2"0 e _2 K 2"0 < r'|p|e C ' |p|2 



(2fc- 1)! 

fe=ln o >0 V ; 
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for some new constants T' , c' . Performing the limit N — > oo in (A1.7), we see that F(ip,p) satisfies 
the same bound 

|F(V>,f>)| < r'|p| e c ' |p|2 , for |Imp| < a, (A1.8) 

for some a small enough, so that F(ip,p) is analytic (in p) in a strip and therefore h(ip, rj) is Borel 
summablc (in rj). 



Appendix A2. Comparison with the method of [GG1] 

In this Appendix we briefly discuss how the function h constructed in [GG1] can be identified with 
the Borel summable function of Proposition 1. By the uniqueness in the class of Borel summable 
functions, it is enough to prove that also the function h of [GG1] is Borel summable. 

We begin by reviewing the differences of the construction envisaged in [GG1] with respect to that 
of [GG2]. Trees, labels and clusters are defined in the same way as in items (a) to (h) of Section 
2. What changes is the definition of the propagators, which is iterative. By writing x — u> ■ v>i, 
Vi 7^ 0, we set gj, ' = l/x 2 , and, for k > 1, 

with (x; if) defined as 

M^{x;rf) =]TVt(x; 77), (A2.2) 

T 

where the sum is restricted to the self-energy clusters T with scale nr > n + 3, where n is such 
that jn-i < M < 7nj an d the self-energy value is given by 

v ^-^=-Jnu{ n ^)( n ^ 

1 v n ' vev(T) ieA(T) 

Note that k labels the iterative step and, in principle, it has no relation with the scale n of 
x. However the matrices (x; if) are obtained from resummations of self-energy clusters with 
height up to k — n (for definitions and details we refer to [GG1], where the self-energy clusters 
are called self-energy graphs). Hence they stop flowing at k — n if x is on scale n: this means that 
M^(x;rj) = M^(x;r]) as soon as k > n. 

Then h' - ' will be expressed in terms of trees as in (2.2), if Y&\(0) is defined as in (2.4), with 
gf"^ replaced with g[ N \ and is obtained as the limit of as N — > 00. 

Let us consider for simplicity the case of polynomial perturbations. Then we can proceed as in 
Sec. 4, and prove by induction the bound 

\\(M^) B (x; P )\\<r\p\e dkM , for alia;, (A2A) 

for all p complex. Supposing inductively the bound (A2.4) we obtain that the Borel transform of 
gf^ can be bounded as 

\\(g [k] )B(x; P )\\ < l 4e c *W. (A2.5) 
x z 

This is trivial for k = 0, as (g^)B(x;p) = p/x 2 , while it follows from (A2.4) for n > 1 by the 
inductive hypothesis (and it can be proved as the analogous bound in Sec. 4). Therefore we can 
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write M^ N \x; if) according to (A2.2), and its Borel transform (M^)b(x;p) can be computed and 
bounded as done in (4.5) and (4.8), so that at the end the bound (A2.4) is obtained for k — N . In 
particular the bounds (A2.5) on the propagators imply that admits the same bound (4.10) 
found in Sec. 4. Therefore we can take the limit for N — > oo, and Borel summability for h follows. 
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